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We elucidate effects of solar density fluctuations on neutrino propagation through the Sun. Using 
data from the recent solar neutrino and KamLAND experiments we provide stringent limits on solar 
density fluctuations. It is shown that the neutrino data constrains solar density fluctuations to be 
less than j3 = 0.05 at the 70 % confidence level, where j3 is the fractional fluctuation around the 
value given by the Standard Solar Model. We find that the best fit to the combined solar neutrino 
and KamLAND data is given by j3 = 0. 

PACS numbers: 14.60.Pq, 26.65. +t, 96.60.Jw, 96.60.Ly 
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I. INTRODUCTION 



The evolution of the Sun is modeled using a standard 
set of input parameters describing the physics of the in- 
terior, including thermodynamic properties (equation of 
state), energy transfer through solar matter (opacity) 
and the rates of the nuclear reactions that power the Sun 
(astrophysical S-factors). The resulting detailed models 
of the Sun predict temperature, density, composition pro- 
files, and neutrino fluxes coming from the nuclear fusion 
reactions. Recently the Standard Solar Model (SSM) has 
been enjoying a tremendous success. Its predictions (see 
QI Q ) have withstood the observational and experimen- 
tal tests very well. For example, SSM predicts frequen- 
cies of the pressure (p-mode) vibrations which can be 
observed on the solar surface. Observation and analy- 
sis of these oscillations, called helioseismology, provide 
detailed information about the Solar interior (for a thor- 
ough introduction to helioseismology and review of the 
observational data see Ref. 03)- The helioseismological 
information can be inverted to obtain the sound-speed 
profile throughout the Sun. The speed of sound in the 
Sun is determined by combining the density and temper- 
ature profiles which arc predicted using the Standard So- 
lar Model. The sound speeds of solar models that include 
element diffusion agree with helioseismological measure- 
ments to better than 0.2 % 0,03. Another test of SSM 
was achieved by the neutral-current measurements at the 
Sudbury Neutrino Observatory (SNO) |||. These mea- 
surements yield a total (all flavors) 8 B solar neutrino flux 
which is in very good agreement with the Standard Solar 
Model predictions. Identifying the precise nuclear reac- 
tions that power the Sun would be another valuable con- 
straint on the SSM It is generally believed that the 
nuclear fusion reactions that power the Sun take place in 
the so-called pp cycle. The recent data also provides a 
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stringent limit (less than 7 %) to the amount of energy 
that the Sun produces via the CNO fusion cycle Q. 

Given these recent successes of the Standard Solar 
Model and the quality of the experimental data currently 
being taken, perhaps the time has come to test some of 
the other aspects of the model. One such test, namely 
looking for fingerprints of the solar density fluctuations in 
solar neutrino spectra was alre ady considered extensively 
starting in the early 1990's 0, E3, E3, E3, E3, E3, EH . 
but at that time the solar neutrino data were not accu- 
rate enough to make a definitive statement. Such fluc- 
tuations may provide an additional probe of the physics 
of the deep solar interior For example, evolution 

theories of stars, in addition to p-modes, also predict the 
existence of buoyancy-driven gravity modes (g- modes). 
Whether g-modes are actually excited in the Sun is an 
open question. These modes are exponentially damped 
in the convective zone of the Sun; unlike the p-modes it 
is not possible to observe the resulting very small am- 
plitudes of the g-modes on the surface of the Sun using 
current techniques even if they are indeed excited in the 
Sun. It is argued that g-modes in the Sun can be ex- 
cited by turbulent stresses in the convective zone or 
by magnetic fields in the radiative zone [l7|. Although 
the latter hypothesis requires large magnetic fields in the 
Sun, a magnetic field as large as 10 7 G in the radia- 
tive zone seems to be permitted by the helioseismic data 
|l9j . It should be emphasized that temperature fluctu- 
ations associated with the large-amplitude g-mode os- 
cillations would have significantly reduced the neutrino 
flux produced in the core of the Sun. Observational up- 
per limits on the surface velocity amplitudes (as well 
the direct measurement of the total solar neutrino flux 
at SNO) rule out such large-amplitude g-mode oscilla- 
tions 201. The scenario we discuss here is the possibility 
of smaller-amplitude g-mode oscillations or some other 
mechanism causing fluctuations in the density profile of 
the Sun affecting neutrinos only through their interaction 
with the matter. (For an alternative mechanism for pro- 
ducing density fluctuations in the Sun through magneto- 
hydrodynamic waves see Ref. [2l|]). At the very least it 
is important to investigate what limits the solar neutrino 
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data (supplemented by the constraints of the KamLAND 
reactor neutrino measurements) would place on the size 
of such fluctuations. 

Our goal in this paper is to revisit the subject of density 
fluctuations and investigate limits placed by the recent 
solar and reactor neutrino experiments on the amount 
of fluctuation. Our preliminary attempts to use earlier 
SNO data were presented in Ref. |22j . In the next section 
we review the formalism and summarize our method for 
doing the global analysis. In Section 3 we present our 
results and discuss their implications. 



II. SOLAR DENSITY FLUCTUATIONS 

We will assume that the electron density N e fluctuates 
around the value, (N e ), predicted by the SSM 

N e (r) = (l+(3F(r))(N e (r)), (1) 

and that the fluctuation F(r) obeys 

(F(r)) = 

(F(n)F(r 2 )) = fa 

(F( ri )F(r 2 )F(r 3 )) = 

(F( ri )F(r 2 )F(r 3 )F( ri )) = (fa fa + f 13 f 2i + / 14 / 23 ) 

; (2) 

where /y = / ( | r j — r^|) gives the correlation between fluc- 
tuations in different places. In Eq. Q we can interpret 
the quantity (3 as the fraction of the fluctuation around 
the density given by the SSM. Throughout the current 
paper, we consider the case of delta-correlated (white) 
noise: 



f(r) = 2r6(r), 



(3) 



with the correlation length r as a parameter. We discuss 
the limitations of this assumption later in this section. 
The Hamiltonian describing neutrino evolution in matter 
can be written as as sum of two terms: 



H = Hq + B(r)M, 



(4) 



where Hq is the standard Mikheev-Smirnov-Wolfenstein 
(MSW) Hamiltonian [23| (for a brief review see 0) with 
the average electron density (-/V e ), B(r) is the fluctuating 
c-number and M is a constant operator. For the two- 
flavor mixing Hq governs the time-evolution in the SSM 
density: 
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Sm 2 f ((t) - cos 29 
AE \ sin 29 



sin 29 
-(C(i) - cos 20) 



where 



C(r) = ^fVe(r)>, 



(5) 



(6) 



6 is the vacuum neutrino mixing angle, Sm 2 is the differ- 
ence in the squared masses of the two neutrino species, E 
is the neutrino energy, and (N e ) is the averaged electron 
number density given by the SSM. In Eq. (jSJ), v x is an 
arbitrary combination of and v T |25l | . The fluctuating 
term B(t) is given by 



B (r) = 0F(r)^(N e (r)}, 



(7) 



and M = C3. Using Eq. J2J) one can show that the 
fluctuation of B(r) satisfies the conditions 



<i*(ri)> 
(B(n)B(r 2 )) 
(B(n)B{r 2 )B(r 3 )) 
(B(ri)B(r2)B(r 3 )B(r 4 )) 



= 



2 f 

a fa 


a (fa fa + fa fa + fa fa) 

(8) 



where 



a(r) = ~^P(Ne 



(9) 



It was shown in Refs. 0, ^3 that with the white noise 
assumption of Eq. @ the fluctuation-averaged neutrino 
density matrix satisfies the equation 
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( P (t)) = ~a 2 r[M, [M, (p(t))}] - i[H (t), ( P (t))}. (10) 



For the two- flavor case, depicted in Eqs. JSJl and J7J), after 
calculating the commutators Eq. H1()J) can be re-written 
as a 3 x 3 matrix equation [!J H3 
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where the individual elements of the density matrix are 

1 



(12) 



In these equations Vf is the probability amplitude for the 
neutrino flavor /, and we introduced the definitions 



A(t) 



AE 



(C(t) -cos 26>), D 



dm 
~4E 



and 



GUN e (r)) 2 p 2 T. 



sin 26>, (13) 



(14) 



We numerically solved Eq. I|ll(l using the technique 
developed in Ref. [2(| which we summarize. We first 
consider the constant density case. If we represent the 
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Electron Neutrino Survival Probability with Fluctuations 
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FIG. 1: Mean survival probabilities for the SSM density pro- 
file and sin ^12 = n/6 calculated as described in the text. 
The correlation length is chosen to be 10 km and the prob- 
abilities are plotted for the percentage fluctuation values of 
(3 = (solid line), (3 = 0.015 (dashed line),/? = 0.03 (dotted 
line),/3 = 0.06 (dot-dashed line). 



column vector in Eq. I|llf) by lZ(t) and the matrix govern- 
ing its time-evolution by K (dIZ/dt = JC1Z), for constant 
density one can write 



K{t) = cxp [K,t]K(t = 0) = U(t)K(t = 0). 



(15) 



The exponential in Eq. I|15f) can be calculated using 
Cayley-Hamilton theorem which states that for any ma- 
trix M, the eigenvalue in the characteristic equation can 
be replaced with the matrix itself. For a 3 x 3 matrix we 
get: 



,M 



= a I + aiM + a 2 M 2 



(16) 



The coefficients can be found by by substituting eigen- 
values Xi of M into 



a + Ajai + X i a 2 . 



(17) 



In a medium with varying density (where the variables 
A and k are changing), time evolution can be calculated 
by dividing neutrino path into small intervals in which 
potential can be approximated as a constant. Then the 
total evolution operator is product of evolution operators 
for all intervals: 



u = u N u N „ 1 u 2 u 1 



(18) 



Care must be taken to ensure that the intervals chosen 
are smaller than the oscillation length. 

To give a feeling about the solutions of Eq. l|llfl we 
show in Figure ^mean survival probabilities obtained by 
solving it numerically as described above for different val- 
ues of P and a correlation length of r = 10 km. From 
this figure one observes that the effect of the fluctuations 



FIG. 2: The effect of source averaging on various components 
of the mean survival probability without (/3 = 0) and with 
(f3 = 0.3) fluctuations. The SSM density profile of Ref. Q, 
the mixing angle of sin 812 = tt/6, and the correlation length 
of 10 km are used. The solid lines are the survival probabili- 
ties without source averaging. The dashed lines represent the 
source- averaged survival probabilities of 8 B neutrinos and the 
dotted-line those of the pp neutrinos. 



is more dominant when the neutrino parameters and the 
average density are such that neutrino evolution in the 
absence of fluctuations is adiabatic. For the two-flavor 
case one can write the electron neutrino survival proba- 
bility as m 113 



- + -cos20cos20 M (l - 2|-0 2 | 2 ) 



(19) 



where \ip 2 \ 2 is the probability of observing the second 
matter eigenstate on the surface of the Sun (calculated 
with the initial condition that the neutrinos start in the 
first matter eigenstate) , also known as the hopping prob- 
ability, and the matter angle is given as 



cos 26 M = —A/y/ A 2 + D 2 . 



(20) 



In Eq. I|19fl cos 29m is the matter angle where the neu- 
trinos are created. Note that the cos2f?M is zero at the 
MSW resonance. If the neutrino propagation is adia- 
batic when j3 is set to zero, the hopping probability is 
zero. However when the fluctuations are turned on they 
cause some hopping, yielding a small but non-zero l^l 2 , 
which grows with (3. Consequently when 7^ 0, the ab- 
solute value of the second term on the right-hand side of 
Eq. (|19|) is always less than its value when (3 = 0. If the 
value of Sm 2 / E is such that neutrinos are produced at 
the MSW resonance density then the cosine of the initial 
matter angle is zero, and Eq. (|T$|) predicts a survival 
probability of 1/2 no matter what the value of (3 is. For 
smaller values of Sm 2 / E the initial value of the matter 
angle is negative, yielding a higher survival probability as 
compared to the (3 = case. For larger values of Sm 2 /E 
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2x2 Solar Neutrino Oscillations with Fluctuations 
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FIG. 3: Source-averaged survival probabilities for 8 B neutri- 
nos with and without solar-density fluctuations detected at 
the location of SNO during the day and the night. The SSM 
density profile of Ref. QJ and the correlation length of f km 
are used. 



the initial value of the matter angle is positive, yielding 
a lower survival probability as compared to the = 
case. This behavior is clearly evident in Figure ^ 

There are two constraints on the value of the corre- 
lation length. In averaging over the fluctuations we as- 
sumed that the correlation function is a delta function 
(cf. Eq. 10). In the Sun it is more physical to imag- 
ine that the correlation function is like a step function of 
size r. Assuming that the logarithmic derivative is small, 
which is accurate for the Sun, delta-correlations are ap- 
proximately the same as step-function correlations if the 
condition 



t < sin 20 



is satisfied llC 



Sm 2 
~2E 

Eq. (|21|l can be rewritten as 
E(MeV) 



r(km) < 3.95 x 10" 



sin 20 Sm 2 {eV' 



(21) 



(22) 



A second constraint on the correlation length is provided 
by the helioseismology. Density fluctuations over scales 
of ~ 1000 km seem to be ruled out HEHH. On the 
other hand current helioseismic observations are rather 
insensitive to density variations on scales close to ~ 100 
km [13. 



III. RESULTS AND DISCUSSION 

In our analysis we used a covariance approach the de- 
tails of which are described in Ref. [2^ . We use 93 data 
points in our analysis; the total rate of the chlorine exper- 
iment (Homestake \3 \) . the average rate of the gallium 
experiments (SAGE[3l|, GALLEX [H, GNO |33|), 44 
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FIG. 4: Allowed regions of the neutrino parameter space 
with solar-density fluctuations when only solar neutrino ex- 
periments (chlorine, all three gallium, SNO and SK experi- 
ments) are included in the analysis. The SSM density profile 
of Ref. Q| and the correlation length of 10 km are used. 
The case with no fluctuations (/3 = 0) are compared with re- 
sults obtained with the indicated fractional fluctuation. The 
shaded area is the 70 % confidence level region. 90 % (solid 
line), 95 % (dashed line), and 99 % (dotted line) confidence 
levels are also shown. 



data points from the SK zenith-angle-spectrum fcyj, 34 
data points from the SNO day-night-spectrum J35j and 
13 data points from the KamLAND spectrum |3(|. We 
took into account the distribution of the neutrino sources 
in the Sun. Specifically we divided the Sun into several 
shells, calculated the survival probability numerically for 
neutrino paths and averaged the survival probabilities 
over the initial source distributions. Similarly we consid- 
ered the effects of the matter density of the Earth (the 
day-night effect) by solving neutrino evolution equations 
numerically. We illustrate the effects of source averag- 
ing (with and without fluctuations) in Figures [21 and [3] 
In Figure |21 we show the source-averaged survival prob- 
ability separately for the pp and 8 B neutrinos. In this 
figure the solid line represents the survival probability 
of the neutrinos coming from a single point (the center 
of the Sun) contrasted to the source-averaged cases. For 
= 0, for a given 5m 2 , neutrino energy, and mixing angle 
neutrino evolution is adiabatic for the region of interest 
shown in the graph. In contrast, for ^ 0, one has a non- 
zero, but small hopping probability. When the neutrinos 
are created over a finite-size region (instead of a single 
point) source-averaged survival probabilities do not coin- 
cide at large 8m 2 /E values with the point-source survival 
probability unlike the = case. Fluctuations also re- 
duce the effect of the averaging. In Figure we present 
the source-averaged survival probabilities detected at the 
location of SNO for 8 B neutrinos with and without solar- 
density fluctuations. To calculate this figure we used the 
day and night live-time information from Ref. |37| . One 
again observes that fluctuations smoothen the survival 
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2x2 Solar Neutrino & KamLAND with Fluctuations 
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2x2 SK + SNO + KamLAND Neutrino Data with Fluctuations 
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FIG. 5: Allowed regions of the neutrino parameter space 
with solar-density fluctuations when the data from the solar 
neutrino experiments (chlorine, all three gallium, SNO and 
SK experiments) and the KamLAND data are included in 
the analysis. The SSM density profile of Ref. Q and the 
correlation length of 10 km are used. The case with no fluc- 
tuations (f3 = 0) are compared with results obtained with the 
indicated fractional fluctuation. The shaded area is the 70 % 
confidence level region. 90 % (solid line), 95 % (dashed line), 
and 99 % (dotted line) confidence levels are also shown. 

probabilities. 

We next turn our attention to parameter-space 
searches. In Ref. [2!| from a global analysis of solar neu- 
trino and KamLAND data we found for electron neutrino 
oscillations into another active flavor, the best fit values 
of tan 2 f?i2 ~ 0.46 for the mixing angle between first and 
second generations, tan 2 #13 ~ for the mixing between 
first and third generations, and Sm^ ~ 7.1 x 10 -5 eV 2 . 
Other groups doing similar analyses found very similar 
best fit values j33- In Ref. J2|| /3 was taken to be zero. 
We find that non-zero values of (3 reduces the size of 
the allowed region as shown in Figure 0] In this fig- 
ure allowed regions of the neutrino parameter space are 
shown for different values of (3 when all the solar neu- 
trino experiments (chlorine, all three gallium, SNO and 
SK experiments) are included in the analysis. Although 
the minimum value of \ 2 is achieved when f3 = 0, one 
observes that for values of (3 as large as 0.07 allowed re- 
gions remain even at the 70 % confidence level. Thus we 
conclude that solar neutrino data alone does not signif- 
icantly constrain the fluctuation parameter. Note that 
for larger values of (3 the region with larger values of 6m 2 
are no longer allowed. Since KamLAND data favor larger 
values of r5m 2 , incorporating KamLAND results dramat- 
ically reduces the allowed region in parameter space as 
shown in Figure Again the he minimum value of \ is 
achieved when (3 = 0. However, in contrast to the calcu- 
lation presented in Figure one can put stringent limits 
on the amount of fluctuation. We find that (3 < 0.05 at 
the 70 % confidence level, (3 < 0.06 at the at the 90 % 
confidence level, and (3 < 0.07 at the at the 95 % confi- 



FIG. 6: Allowed regions of the neutrino parameter space 
with solar-density fluctuations when the data from the solar 
neutrino experiments detecting only the high energy 8 B neu- 
trinos (SNO and SK experiments) and the KamLAND data 
are included in the analysis. The SSM density profile of Ref. 
PU and the correlation length of 10 km are used. The case 
with no fluctuations (/? = 0) are compared with results ob- 
tained with the indicated fractional fluctuation. The shaded 
area is the 70 % confidence level region. 90 % (solid line), 95 
% (dashed line), and 99 % (dotted line) confidence levels are 
also shown. 

dence level. 

Let us recall that both calculations are done for a value 
of r = 10 km for the correlation length. Note that only 
the combination /3 2 r enters into the calculation (cf. Eq. 
1)14(1. Hence the (3 values can be scaled by adjusting the 
value of r and for larger correlation lengths the limits 
quoted above get tighter. However one cannot consider 
arbitrarily large values of the correlation length. Our for- 
mulation of the problem (the delta-function correlation 
approximation) becomes unrealistic for larger values of 
t as we illustrated in Eq. 1(22(1 . In we insert the best 
fit (minimum \ 2 ) values of 8m 2 ~ 7.1 x 10~ 5 cV 2 and 
tan 2 6 ~ 0.46 into Eq. (55) we find 

r(km) < 6£(MeV). (23) 

Hence for lower-energy (pp) neutrinos the reliable cor- 
relation lengths are smaller then 10 km. However for 
higher energy ( 8 B) neutrinos one can safely consider 
longer correlation lengths. For both SK and SNO the 
energy threshold is ~ 5 MeV for which Eq. 1(23(1 yields 
t < 30 km. To explore this feature we repeat our anal- 
ysis considering only SK and SNO data together with 
the KamLAND data. We show the allowed regions of 
the parameter space in Figure El In this figure for better 
comparison to Figure we took r to be 10 km, however 
even r = 20 km would be reasonable. The best fit is still 
with = 0. We find that (3 < 0.07 at the 70 % confi- 
dence level with r = 10 km. This limit would scale down 
to 13 < 0.07/V2 ~ 0.05 at the 70 % confidence level with 
r = 20 km. 
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Effect of KamLAND on % 
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FIG. 7: Projection of the global A\ function on the frac- 
tional fluctuation parameter /3. The solid line represents the 
calculation using only the solar neutrino data. The dotted 
line is the calculation when the KamLAND data are included 
along with the solar neutrino data. For comparison we also 
show the calculation with solar-neutrino data sensitive only 
to the 8 B neutrinos (SK and SNO) combined with the Kam- 
LAND data. 



role to constrain (3. As the KamLAND statistics improve 
in the near future we expect to improve our limits on the 
fractional fluctuation. 

In this paper we ignored fluctuations of another kind, 
namely magnetic field fluctuations impacting neutrino 
evolution. It has been shown that if neutrinos have siz- 
able magnetic moments they can interact with the trans- 
verse magnetic fields [3|j undergoing a spin-flavor pre- 
cession [4(|. If the magnetic field is n oisy the spin-flavor 
precession will be impacted |^ ITial4lll42 |. However a rel- 
atively conservative assumption about the maximal size 
of the solar magnetic field places the spin-flavor reso- 
nance at values of 8m\ 2 0] which is at the region of the 
neutrino parameter space ruled out by the KamLAND 
experiment. (For an alternative approach see Ref. [i3|L 
Hence one can conclude that solar magnetic field fluctu- 
ations do not play a role in the solar neutrino physics 
through the spin-flavor precession mechanism. 

We would like to point out that besides in the Sun 
(and other stars) neutrinos interact with dense matter in 
several other sites such as the early universe, supernovae, 
and newly-born neutron stars and neutrino interactions 
with a stochastic background may play an even more 
interesting role in those sites Along those lines a 

preliminary analysis of the effects of random density fluc- 
tuations on matter enhanced neutrino flavor transitions 
in core-collapse supernovae and implications of such fluc- 
tuations for supernova dynamics and nucleosynthesis was 
given in Ref. jlq . 
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